The role of non-local exchange in the electronic structure of correlated oxides 
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We present a systematic study of the electronic structure of several prototypical correlated 
transition-metal oxides: VO2, V2O3, Ti203, LaTiOa, and YTiOa. In all these materials, in the 
low-temperature insulating phases the local and semilocal density approximations (LDA and GGA) 
of density-functional theory yield a metallic Kohn-Sham band structure. Here we show that, without 
invoking strong-correlation effects, the role of non-local exchange is essential to cure the LDA/GGA 
derealization error and provide a band-structure description of the electronic properties in quali- 
tative agreement with the experimental photoemission results. To this end, we make use of hybrid 
functionals that mix a portion of non-local Fock exchange with the local LDA exchange-correlation 
potential. Finally, we discuss the advantages and the shortcomings of using hybrid functionals for 
correlated transition-metal oxides. 

PACS numbers: 71.20.-b,71.30.+h,71.15.-m 



I. INTRODUCTION 

Nowadays, the standard model of electronic struc- 
ture calculations is based on density-functional theory 
(DFT) in the Kohn-Sham (KS) formalism!. The DFT- 
KS scheme, also in its simplest approximations like the 
local-density approximation (LDA) 1 or the generalized- 
gradient approximation (GGA)2, is generally highly suc- 
cessful in a very large variety of applications. Thus, when 
a LDA (or GGA) KS band structure turns out to be 
metallic in an insulating compound (e.g. in a transition- 
metal oxide), the result is often interpreted as a direct 
indication for strong electron correlation effects in the 
material, and a failure of the band-structure picture. A 
possible strategy that has been followed in the literature 
to overcome these difficulties is to resort to model ap- 
proaches, like the multiband Hubbard model. In the 
LDA+U approach^, LDA band structures are supple- 
mented by an on-site Coulomb interaction (the Hubbard 
U) acting only on the "correlated" subset of the electronic 
degrees of freedom. In a higher level of theory, dynamical 
mean- field theory (DMFT)i, the Hubbard model is fur- 
ther mapped onto an Anderson impurity model, which 
can be then solved with different techniques, allowing for 
the description of dynamical effects beyond LDA+U. In 
those cases LDA is claimed to be inadequate to capture 
the strong interactions taking place between correlated 
electrons in partially filled d (or /) shells, which give rise 
to narrow bands in the solid. 

However, we remark that DFT is a ground-state the- 
ory and KS band structures are not meant to describe 
the electronic excitations measured in photoemission, 
which also define the fundamental band gap of an in- 
sulator. Moreover, LDA lacks the derivative discontinu- 
ity of the local KS exchange-correlation (xc) potential 5 
V xc (r) and suffers from a severe delocalization error—, 



which is particularly relevant for localised d and / elec- 
trons. The underestimation of the fundamental band 
gaps in sp semiconductors is well knowni£, and under- 
stood in terms of self-energy corrections at the GW level 
of approximation^^. This underestimation sometimes 
may lead to metallic band structures also in "weakly cor- 
related" small-gap semiconductors, like in germanium. 

Here we consider several prototypical correlated 
transition-metal oxides that have been studied by other 
methods in the past years: VO2, V2O3, Ti203, LaTi03, 
and YTi03. In all the low-temperature insulating phases 
of these materials, KS-LDA yields a metallic band struc- 
ture. We thus address the following question: Is this just 
a result of the inadequacy of LDA to deal with strong 
correlations or, rather, is this finding related to the sys- 
tematic KS-LDA underestimation of band gaps that oc- 
curs also in "weakly correlated" semiconductors and get 
enhanced in the oxides? To answer this question we 
make use of a generalized Kohn-Sham (gKS) scheme^, 
where the local KS xc potential V xc (r) is replaced by 
a spatially non-local V xc (r,r'). In these hybrid func- 
tionals, the non-local Fock exchange potential is mixed 
with the local (LDA or GGA) KS xc potentially. In 
these approaches, the difference with LDA results stems 
solely from the non-local exchange term. This is not 
supposed to improve the description of electronic corre- 
lation. Therefore, if strong correlations are responsible 
for invalidating the KS-LDA description of these insula- 
tors, then one should expect to find the same problems 
in the gKS scheme. We will instead show how, without 
invoking strong-correlation effects, the role of non-local 
exchange is essential to cure most of the LDA delocaliza- 
tion error and provide a band-structure description of the 
electronic properties of several transition-metal oxides in 
qualitative agreement with the experimental results. 

The paper is organized as follows. In Sec. [H]we briefly 



2 



introduce the generalized Kohn-Sham (gKS) scheme and 
the hybrid-functional paramctrization that we use, and 
compare them with standard methods that treat elec- 
tronic correlations. In Sec. IHII we present and discuss the 
results that we have obtained for the various transition- 
metal oxides, also across their metal-insulator phase tran- 
sitions (MIT)i£. Finally, in Sec. HVI we draw our conclu- 
sions on the basis of these results and discuss the advan- 
tages and shortcomings of the use of hybrid functionals in 
the description of the electronic properties of correlated 
oxides. 



II. METHOD 

In the gKS scheme the xc potential is generalized to 
be non-local in spaced, contrary to ordinary KS where 
it is local. In this framework, a common choice is to 
write the gKS non-local potential V xc (r,r') as a sum of 
a non-local exchange term and a local KS correlation 
potential^: V xc (r,r') = V x (r, r') + V c (r). For V c (r) a 
LDA or GGA expression is adopted, while the non-local 
exchange term V x (r,r') is obtained mixing the KS local 
exchange V x (r) with a fraction a of the non-local Fock 
operator V F (r, r') — —j(r,r')v(r,r'), built using gKS 
orbitals in the one-particle density matrix 7 [v(r,r') — 
l/\r — r' I is the Coulomb interaction]: 

V x (r,r') = aV x F (r,r') + (1 - a)V x (r). (1) 

This construction is often justified from considerations 
based on the adiabatic-connection formula, which fix the 
mixing parameter to be a = 0.25. 1 — The Coulomb in- 
teraction can be then split into a sum of a long-range 
and a short-range term v(r, r') — v LR (r, r') + v (r, r'), 
respectively as: 

t,(ry) = erf ^ |r "^ l) + 1 ~f (M|r ~ r,|) . (2) 
\r — r I \r — r \ 

Here both the choice of the erf error function and the 
parameter (i are arbitrary. Replacing the bare interac- 
tion v by either term, the same separation is obtained 
in the local and non-local exchange potentials appear- 
ing in Q: V x F (r,r') = V FLR (r,r') + V F ' SR (r,r') and 
V x (r) = V£ R {r) + V^f H (r). Assuming that the effects 
of V F,LR (r,r') and V x LR (r) compensate each other, the 
following approximation for (fTJ) is introduced: 

V x (r,r<) = aV x F ' SR (r,r') + (l-a)V x SR ir) + V x LR (r). (3) 

Regrouping the various terms contributing to the xc gKS 
potential, one finally finds: 

V xc (r, r') = a[V F ' SR {r, r') - V x SR (r)] + V xc (r), (4) 

where the correction to the KS local potential V xc (r) 
(with the original Coulomb interaction) stems entirely 
from the first term in the r.h.s. The final approximation 
depends on the two parameters a and fi. With fi = 0, 



a = 0.25, and using the PBE GGA xc potential for the 
local part V xc (r), one finds the PBEO approximation. 13 
In the HSE06 hybrid functional^, instead, the value of 
fx = 0.2 A is obtained by numerically fitting the results 
against a benchmark set of data. 

More in general, both the a and fi parameters play the 
role of effective screening of the Coulomb interaction^. 
The non-local potential (fTJ) with fj, = can be alterna- 
tively understood as a static approximation to the many- 
body GW self-energy 9 -. By identifying 1/a as an effec- 
tive static dielectric constant e, aV F can be seen as a 
screened exchange potential, while the local part of (TTJ) 
acts as an approximation to the Coulomb hole term&H 
In fact, varying a between and 1 in Eq. (TTJ), one 
in practice interpolates between the KS underestimation 
and the Hartree-Fock (HF) overestimation of band gaps, 
with the possibility to get close to the experimental re- 
sults. Moreover, the use of a finite value for fi, together 
with neglecting the corresponding long-range terms in 
Q , efficiently acts as a further screening of the Coulomb 
interaction . 14 i 18 By increasing the value of /i, the two- 
point distance \r — r' | beyond which the Coulomb in- 
teraction is cut off becomes shorter. In fact, screening 
the long-range Coulomb interaction is crucial for many 
properties in bulk materials 19 , especially for small-gap 
and metallic systems. While for /1 = 00 Eq. ^ reduces 
to the KS local potential V xc (r), a finite fi tunes the 
correction to KS stemming from the difference between 
non-local Fock and local exchange terms in . For fixed 
a, increasing fi gives less weight to this correction. 

While LDA suffers from a derealization error, HF in- 
stead is affected by an excess of localization 6 . Hence, the 
inclusion in the functional of a partial contribution of 
non-local Fock exchange leads to a localization correction 
with respect to LDA. This affects in particular d and / 
states, which are more localised and generally suffer from 
a self-interaction error more than s and p states^ Even 
though it is derived from a very different point of view, 
this effect of localization of d and / orbitals is shared 
also by the LDA+U approach. As we will show in the 
following, curing the derealization error of LDA is the 
key to get results in better qualitative agreement with 
experiment. On the other hand, both LDA+U (as well 
as LDA+DMFT) and the hybrid functional (@} depend 
on parameters (like the Hubbard U and J in the former 
and a and fi in the latter), which limits their predictivity 
power. However, in the hybrid functional (j4|) there is no 
need of an additional term for correcting the (unknown) 
spurious double counting, while it is case for LDA+U or 
LDA+DMFT. Moreover, in the hybrid functional all the 
electrons are treated on equal footing. 

Modeling efficiently the static screening of the 
Coulomb interaction without introducing adjustable pa- 
rameters would help to improve greatly the gKS scheme. 
At the same time, this should keep its computational cost 
cheaper than the more sophisticated many-body GW ap- 
proximation, where the dynamical screening is explicitly 
calculated (in the random-phase approximation). In any 
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case, being a static approximation to the many-body self- 
energy, the V xc (r,r') functional in (fj| cannot account for 
dynamical correlation effects, which for instance give rise 
to satellites in the photoemission spectral, and are in- 
stead accounted for, though in different manners, by both 
the GW approximation and the LDA+DMFT approach. 

In the following we will use the hybrid functional form 
(fJJ , as implemented in the VASP computer cod o 22 i 23 . For 
the local V xc {r) part we adopt a LDA xc potential. In the 
case of VO2 we will discuss the effect of different choices 
of the two parameters a and /1, while for the rest of the 
paper we will fix a = 0.25 and /1 = 0.2 A as in the 
HSE06 functional. 



III. RESULTS AND DISCUSSION 



vo. 



VO2 undergoes a twofold phase transition at 340 K21. 
The MIT is accompanied by a lowering of the symmetry 
of the crystal structure, from rutile to monoclinic, with 
a doubling of the unit cell and a dimerization of V atoms 
along the rutile c axis. It has been long debate d 25 i 26 
which of the two aspects, the electronic or the structural 
change, is the key to drive the phase transition. 

In VO2 both LDA 2 - and standard single-site 
DMF T 28 i 29 are unable to get the insulating band gap, 
while LDA+U 28 i 30 has problems with the metallic phase 
and gives an ordered magnetic phase for the insulator, 
contrary to the experiment. The deficiencies of single- 
site DMFT have been corrected by its extension to clus- 
ter DMF T 31 i 32 , where the local impurity is taken to be a 
V dimer instead of a V atom, as in single-site DMFT. On 
the other hand, parameter-free GW calculations^ 3 - have 
shown that in the insulating phase KS-LDA wavefunc- 
tions are not a sufficiently good approximation to quasi- 
particle (QP) wavefunctions. The LDA error is due to an 
excessive derealization of the KS wavefunctions at the 
Fermi level. They turn out to be too isotropic, underesti- 
mating the effect of the V dimerization along the c axis 2 ^, 
and the corresponding bonding- antibonding splitting of 
the V dig states. Once this LDA error is corrected by us- 
ing better QP wavefunctions, as obtained in a restricted 
self-consistent GW scheme, the results correctly repro- 
duce the electronic properties of both phases 3 - 3 - - — and 
also show that the satellite in photoemission spectrum 
of the metallic phase is related to a neutral (plasmon) 
excitation visible in the loss function 3 -^. 
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FIG. 1. (Color online) Densities of states of insulating VO2 
obtained with the hybrid functional of Eq. ((4| according to 
different choices of the mixing parameter a and the screening 
parameter fi. In (a) the topmost panel corresponds to LDA. 
In the other panels /i = 0.2 A -1 and a is raised up to 0.5. In 
(b) a is fixed to 0.25 and the results are obtained for different 
values of \i. In the topmost panel the PBE0 result is retrieved. 
Here, and in all the following figures, the Fermi energy for 
insulators is set in the midpoint of the band gap. 



A very recent calculation^ within the HSE06 hybrid 
functional obtained a gap in density of states (DOS) of 
the insulator, also concluding in favor of the structural 
distortion as the key to explain the insulating gap. Here, 
more in detail, we start our investigation by analysing 
the performance of the hybrid functional Q according to 



different choices of the values of the mixing parameter a 
and the screening parameter /i. We use the experimental 
crystal structures for the two phase s 38 ' 39 , a 6x6x6 k- 
point grid for the insulator and a 12x12x12 one for the 
metal. As we can see in Fig. [TJ the hybrid functional (|4]) 
correctly yields a gap in the insulating phase of VO2 for 
many choices of the two parameters, correcting the LDA 
error [which is retrieved for a = 0, see topmost panel 
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TABLE I. Values of the fundamental band gap, the p band- 
width, the Op - V Oi 9 separation, and the V a lg bandwidth 
in insulating VO2 depending on the different choices for the 
a and fi parameters in the hybrid functional of Eq. @. 



a 




Band gap 


O p 


Op-Vaig 


Vai 9 




(A- 1 ) 


(eV) 


(eV) 


(eV) 


(eV) 







0.00 


6.23 


0.96 




0.125 


0.2 


0.00 


6.54 


1.02 


0.60 


0.250 


0.2 


1.13 


6.85 


0.95 


0.63 


0.500 


0.2 


1.17 


7.58 


0.09 


1.16 


0.250 


0.0 


1.88 


6.86 


0.94 


0.64 


0.250 


0.4 


0.65 


6.74 


0.96 


0.60 


0.250 


0.6 


0.00 


6.62 


0.97 


0.59 



in Fig. [Ha)]. On the other hand, the DOS is highly 
sensitive to the values of a and /1, both for the size of 
the band gap and for the separation between the highest 
occupied V a\ g state and the rest of the O p states at 
higher binding energies (see Tab. IJ). Either reducing a 
[Fig. QJa)] or increasing /i [Fig. QJb)], the gap decreases 
until it disappears for a < 0.125 or /1 > 0.6 A -1 . Keep- 
ing a fixed, the reduction of /i from 0.6 to A -1 has 
mainly the effect of a rigid expansion of the DOS [see 
Fig. Efb)]. On the contrary, fixing fi and changing a 
leads also to larger modifications in the shape of all the 
structures appearing in the DOS, both for occupied and 
the unoccupied states [see Fig. [TJa)]. Thus, the addition 
of non-local exchange to the KS-LDA functional leads 
to an improved qualitative agreement with experiment, 
i.e. a sizable band gap is obtained, without the need to 
invoke strong correlation effects. However, the compar- 
ison with photoemission spectra depends quantitatively 
on the choice that one makes for the values of the two 
parameters a and /i. In fact, they both physically act 
as screening of the Coulomb interaction. However, their 
fine tuning generally depends on an adequate microscopic 
description of the screening in the actual material. 

In Fig. [5] we compare the calculates DOS for both 
phases of VO2 with the experimental photoemission 
spectra—. When comparing with experiments, it is es- 
sential to have bulk-sensitive photoemission data, which 
can be obtained by using high-energy photons as in a 
X ray photoemission (XPS) or, even better, a hard X 
ray photoemission (HAXPES) setup—. In fact, at low 
photon energies, photoemission spectroscopy is mainly 
surface sensitive. However, in these materials the elec- 
tronic properties of the surfaces are generally different 
from the bulk. Here, and in the rest of paper, we use 
\i = 0.2 A -1 and a = 0.25, in agreement with the HSE06 
parametrization. With this choice we find a gap of 1.13 
eV, in excellent agreement with the value of 1.1 eV from 
Ref.— , where a GGA instead of LDA KS local functional 
has been used though. This shows that using either LDA 
or GGA KS functional does not have an influence on the 
DOS. Both results overestimate the experimental band 
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FIG. 2. (Color online) Comparison of the calculated density 
of states of insulating monoclinic (upper panel) and metallic 
rutile (bottom panel) phases of VO2 with the experimental 
photoemission spectra from Ref—. The MIT is correctly re- 
produced by the HSE06 hybrid functional, while the LDA 
DOS is always metallic. 



gap, which is 0.6 eV in the insulator 40 . Nevertheless, the 
HSE06 hybrid functional is able to describe correctly the 
MIT, contrary to LDA, LDA+U, and LDA+DMFT. 



B. V2O3 

Like VO2, V2O3 is a time-honored correlated ma- 
terial, whose great interest is due to its temperature- 
induced MTTi^. At T > 154 K it is a paramagnetic 
metal, while at low temperature it becomes an anti- 
ferromagnetic insulator and undergoes a crystal distor- 
tion from a corundum to a monoclinic structure. The 
phase diagram is made more complex by doping with Cr, 
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which induces a different isostructural MIT to a param- 
agnetic phase. The antifcrromagnetic phase has been 
studied in LDA+U 4 ^ (also followed by a perturbative 
GW calculation 44 ) and LDA+DMFTi 5 ., while the spec- 
tral properties of the metallic phase have been extensively 
analysed in LDA+DMFTi 5 .^ and in GW^. 

Here we take the experimental lattice parameters of 
the pure compound s 50 ! 51 . We consider different mag- 
netic configurations for the insulating phase, in order to 
analyse their influence on the electronic structure and 
compare hybrid functional results with those found in 
LDA+U 43 . In the calculation we used a 6x6x6 grid of k 
points for the insulator, which becomes 10 x 10 x 10 for the 
metal. In the experimental magnetic structure (AFI1) 52 , 
each V atom has one spin-parallel neighbor and two spin- 
antiparallel neighbors in the (distorted) hexagonal planes 
of the monoclinic crystal, while the coupling between 
neighbors in different planes is ferromagnetic [see inset 
to Fig. [3][a)]. We also considered a ferromagnetic (FM) 
order and an alternative antiferromagnetic (AFI2) order, 
in which, with respect to experiment, the interplane mag- 
netic coupling is inverted, while is unchanged inside the 
planes [see inset to Fig. [3jL>)]. LDA always yields a 
metal, regardless of the magnetic structures. Instead, 
HSE06 gives an insulator for all the magnetic configura- 
tions considered. The band gap is 1.80 eV in AFI1, much 
larger than the 0.66 eV experimental optical gap 5 -^. The 
two antiferromagnetic DOS turn out to be very similar, 
while for the FM the DOS [see Fig. [3(c)] is quite differ- 
ent and the band gap reduces to 0.7 eV. In the LDA+U 
calculation 4 -, the FM DOS is half-metallic, and the ex- 
perimental antiferromagnetic structure has a 0.7 eV gap. 
HSE06 consistently overestimates the band gap in all 
the magnetic structures. This result is connected to an 
overestimation of the local magnetic moment, which in 
HSE06 is 1.8 hb for the AFI structures and 1.9 [Ib for 
the FM, whereas experimentally it is 1.2 /is per V atom. 
Similarly to LDA+U, the ground-state total energy dif- 
ference between different magnetic structures is rather 
small, with the experimental AFI1 structure being the 
one with lowest total energy. 

In Fig. 2] we compare the calculated DOS with the 
experimental HAXPES results from Ref^. As in VO2 
(see Sec. HITA"]) . HSE06 is able to reproduce the MIT, 
contrary to LDA. At the same time, it also corrects the 
LDA underestimation of the binding energy of the O p 
states in both phases. This is a clear illustration of the 
advantage of treating all the electrons on equal footing 
(in LDA+U or LDA+DMFT, instead, the position of O 
p states is essentially the same as in LDA). However, 
HSE06 overestimates the band width of the top valence 
dig states, and the band gap in the insulator. 



C. Ti 2 3 

Ti203 has a corundum crystal structure as metallic 
V2O3 and a formal d 1 configuration as VO2. Below 
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FIG. 3. (Color online) Density of states of monoclinic insu- 
lating V2O3, calculated both in LDA and HSE06, according 
to the different considered magnetic structures. The mag- 
netic order is visualized in the insets to each panel, where 
the light blue circles schematically represent V atoms in the 
(distorted) hexagonal planes that characterise the monoclinic 
crystal structure and the arrows display the direction of the 
local magnetic moments. 
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FIG. 4. (Color online) Comparison between the DOS calcu- 
lated in LDA and HSE06 and the experimental photoemission 
spectra from Ref.— . (a) Antiferromagnetic insulating phase 
and (b) paramagnetic metallic phases of V2O3. 
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FIG. 5. (Color online) (a) LDA and HSE06 densities of states 
(upper and bottom panels, respectively), together with the 
DOS projected onto Ti s, Ti d and O p states, for insulating 
Ti203 at room temperature, (b) Comparison of the total DOS 
in LDA and HSE06 with the XPS spectrum from Ref.^i. 



400 K it is insulating and above 600 K it is metallic, 
undergoing a very broad temperature-induced MI T 24 ' 55 . 
Moreover, neither phase is magnetically ordered 5 -^, and, 
contrary to both V0 2 and V 2 3 , the MIT in Ti 2 3 is 
isostructural^,. In fact, by raising the temperature only 
an increase of the c/a ratio is observed, accompanied by 
an increase of the Ti dimer distance along the c axis. 
Here we consider two crystal structures with lattice pa- 
rameters measured at 296 K and 868 K— , where Ti 2 C>3 is 
insulating and metallic, respectively. In our calculations 
we used a 6x6x6 k-point grid for the insulating phase 
and a 8x8x8 one for the metal. 

The presence of Ti dimers along the c axis leads to 
a large bonding-antibonding splitting of the ai g states, 
as in the insulating V0 2 case. Here, without additional 
structural changes, we can directly relate the shortening 
of Ti dimer distance with the increase of the a>i g bonding- 
antibonding splitting, and hence the band-gap opening 
between the bonding a\ g and the e g states in t 2g subband 



of Ti 3d orbitals. Also in the present case, this effect is 
masked in LDA by its underestimation of the anisotropy 
introduced by the Ti dimers. For this reason, the LDA 
DOS remains metallic also at low temperature s 58 ' 59 . The 
relevance of the Ti dimers is underlined also by the fact 
that, as in V0 2 , single-site DMFT is not able to obtain 
the insulating phase, which is instead reproduced by a 
cluster DMFT calculation^ for which the local impurity 
is given by the pair of Ti atoms. 

The effect of the HSE06 corrections over LDA [see Figs. 
03a) and [5] is visible for both phases: i) at the Fermi en- 
ergy, where the top valence band is split off from the bot- 
tom of the conduction states, leading to the gap opening 
in the insulator and to a strong reduction of the spec- 
tral weight at the Fermi energy in the metal; ii) in the 
transfer of spectral weight in the unoccupied Ti i 2g band 
towards the high-energy part of the band; iii) in the in- 
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FIG. 6. (Color online) Comparison between LDA and HSE06 
DOS and Ti s, Ti d and O p projected DOS, calculated for 
metallic Ti 2 3 at T = 868 K. 



creasing of the separation between Ti d states and O p 
states, as a consequence of a rigid shift of the latter. In 
general, HSE06 results do not change the hybridization 
character of the bands with respect to LDA. 

At room temperature HSE06 gives a 0.57 eV gap, 
which is much larger than the 0.11 eV estimate from con- 
ductivity and thermoelectric coefficient measurements 62 . 
However, contrary to HF calculations^ 3 -, spin-polarized 
HSE06 correctly yields a non-magnetic solution also in 
the insulating phase. This result is important since it 
shows that within the hybrid functional the existence of 
a gap is not necessarily linked to the presence of mag- 
netic order, as it often happens in LDA+U. Moreover, 
HSE06 correctly describes also the metallic phase and 
the MIT (see Fig IB]), which is often problematic within 
LDA+U. Thus, these results appear to be two advan- 
tages of HSE06 with respect to LDA+U in non-magnetic 
insulators and metals. 

In the XPS spectrum [at Yiv — 1487 eV, reproduced in 
Fig. [SJb)] 61 for the insulating phase at room tempera- 
ture a satellite is clearly visible at 2.4 eV, between the Ti 
3c? peak at 0.7 eV and the broad O 2p band at 4-10 eV. 
This satellite cannot be obtained with the static hybrid 
functional employed here, but is absent also in the clus- 
ter DMFT calculation^. The HAXPES spectrum (hu = 
5931 eV)£i, which is dominated by the cation s contribu- 
tion as in the vanadium oxide s 36 ' 49 (see Fig. 2]), confirms 
that the satellite present in XPS spectrum is a genuine 
bulk feature of insulating Ti2 03 (instead, experimental 
photoemission results for the metallic phase are not avail- 
able). Overall, the HSE06 results [see Fig. EDJb)] compare 
much better with the experimental spectra than LDA. 



D. LaTi0 3 and YTi0 3 

In a seminal paper, Fujimori et a/. 64 considered a series 
of d 1 transition metal oxides. On the basis of a single- 
band Hubbard model, they explained the opening of the 
band gap, going from metallic VO2 and SrV03 to insu- 
lating YTi03 and LaTi03, as the progressive increase of 
the ratio between the Hubbard U and the d bandwidth, 
which is accompanied by the transfer of spectral weight 
from the quasiparticle peak at the Fermi energy in the 
metals to the Hubbard bands in the insulators. Within 
this view, both LaTi03 and YTi03 are Mott insulators, 
with a gap opening between the lower and the upper Hub- 
bard bands. Thus a band-structure description would 
not be able to reproduce these incoherent atomic-like ex- 
citations in the spectra, and hence the gap. In fact, LDA 
KS band structures are metallic for both compound o 65 ' 66 , 
while LDA+U 6 ^ 6 ! an d LDA+DMFT§§ri2, calculations 
correctly yield a gap. 

At low temperatures both compounds order magnet- 
ically. Below 148 K LaTi03 displays a G-type antifer- 
romagnetic order with a local magnetic moment of 0.57 
Het^. YTi03 is ferromagnetic below 29 K, where the 
magnetic moment is 0.8 (Xb per Ti atomic. It has been 
much debated whether the magnetic properties of these 
compounds can be explained in terms of the formation 
of an orbital liquid or orbital ordering 7 - 3 - - — . 

Here we use the low-temperature experimental crystal 
structures from Refsi 71 ' 76 , and a 4x4x4 grid of k points. 
While LaTiOs is non-magnetic in LDA, in HSE06 the 
local magnetic moment is 0.76 /is, overestimating the 
experimental value of 0.57 [Ib- Also in ferromagnetic 
YTi03, in HSE06 it increases up to 0.84 fis from 0.7 \ib 
in LDA, reaching a similar value than in GGA+U 6 - 7 -. 

Also in these compounds, the HSE06 gives insulating 
densities of states [see Figs. [7[a) andlSJa) for LaTiOs and 
YTi03, respectively]. This is the result of the splitting 
off of a Ti t2 g band from the states crossing the Fermi 
level in LDA. In the case of YTi03 only spin up states 
contribute to the occupied t2 g band. Thus, these results 
for both compounds seem to be in contrast with the tra- 
ditional interpretation of the topmost occupied state as 
an incoherent lower Hubbard band^ 4 -. In fact, the peak 
in the experimental spectra^ 7 - [see Figs. Efb) and [5[b)] 
is matched, at least partially, by this (coherent) Ti ti g 
band. Moreover, within the hybrid functional scheme, 
the opening of the gap with respect to metallic LDA DOS 
is due to non-local exchange. However, contrary to the 
experiment, in HSE06 the band width of these valence 
Ti t2 g states is larger in YTiOs than in LaTiOs. At the 
same time, also the band gap is larger in LaTiOs (1.74 
eV) than in YTi03 (1.41 eV), while the experimental op- 
tical gap is 0.2 eV for LaTiOgS and 0.7 eV for YTiOjZ 9 -. 
In both compounds the band gap is overestimated, while 
the occupied Ti ti g band width is underestimated, also 
for possible dynamical effects 6 ^ that are missing here. 

While, as in all the other compounds, the hybridiza- 
tion between O p states and Ti d states does not change 
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FIG. 7. (Color online) (a) DOS and projected DOS for 
LaTiOa calculated in LDA and HSE06 and (b) comparison 
with experimental photoemission spectra from Ref.—. 
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FIG. 8. (Color online) (a) DOS and projected DOS for YTiOa 
calculated in LDA and HSE06 and (b) comparison with ex- 
perimental photoemission spectra from Ref.—. 



much between LDA and HSE06 [see Figs. [Tfa) and[8Ja)], 
a larger effect is seen here for the unoccupied La / states 
in LaTiOa. In LDA they are located in the middle of the 
conduction band, while in HSE06 they are shifted to the 
upper end of the band. The band gap opening between 
Ti d states and the upshift of the La / states can been 
obtained in the GGA+U approach^ only with the simul- 
taneous use of two (different) Hubbard U values applied 
to Ti d and La / states. Within hybrid functionals this 
result emerges naturally as a consequence of the local- 
ization of these states (the non-local exchange corrects 
the LDA derealization error treating all the electrons 
on equal footing). 

Both LaTiC>3 and YTiC>3 remain insulating also above 
the corresponding (Neel or Curie) temperatures where 
they loose their magnetic order. A spin-unpolarized 
HSE06 calculation would not be able to obtain a gap 
in this case. In agreement with Mott picture^, while the 



long-range magnetic order is not essential to have an insu- 
lator, the electronic spins do matter. In fact, LaTiOa and 
YTiC>3 in the disordered phases are both paramagnetic. 
Similarly, above the Neel temperatures transition-metal 
monoxides are also paramagnetic insulators. For these 
compounds it was recently shown that a calculation, 
based on the self-interaction-corrected (SIC) functional 
and taking explicitly into account the disordered local 
moments of the paramagnetic phase, was able to cor- 
rectly describe the insulating phases^. Thus, a similar 
calculation would be suitable in the (disordered) param- 
agnetic phase also for the present perovskite compounds 
(for which a modified SIC implementation has been re- 
cently used for the magnetically ordered phases^). 
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IV. CONCLUSIONS 

The hybrid fvmctionals employed in the present work 
are not explicitly designed to treat electronic correla- 
tions. Nevertheless, we have shown that the inclusion 
of non-local Fock exchange is essential to cure the fun- 
damental problem of getting metallic band structures in 
Kohn-Sham LDA in the insulating phases of several cor- 
related transition-metal oxides, as those that have been 
discussed here: V0 2 , V 2 3 , Ti 2 3 , LaTi0 3 , and YTi0 3 . 
Analogous results have been obtained for instance also by 
Rodl et al£^ in the series of transition-metal monoxides, 
where hybrid functionals were also used as an improved 
starting point for one-shot GW calculations^. 

Thus, a common conclusion emerges from the study of 
all these correlated transition-metal oxides using HSE06. 
With respect to LDA, besides providing a finite band gap, 
they also correct the position of the O p states. This 
is a clear advantage with respect to other approaches 
(LDA+U, LDA+DMFT, etc.) stemming from treating 
all the electrons on the same footing. Moreover, an- 
other advantage with respect to LDA+U is the consis- 
tent treatment of insulators and metals and the fact 
that hybrid functionals do not yield an insulator together 
with magnetic long-range order, as it is often the case in 
LDA+U. However, the HSE06 parametrization overesti- 



mates the band gap in all the compounds considered here. 
And, in general, the results depend on the choice of the 
parameters (as in LDA+U and LDA+DMFT) used to 
build the functional. Moreover, hybrid functionals miss 
completely dynamical correlation effects, which are es- 
sential for the description of satellites in photoemission 
spectral. Therefore, it is evident that they cannot be 
(and they are not meant to be) the final answer for the 
description of spectral properties of correlated transition- 
metal oxides. However, also in these compounds, they 
demonstrate to be very useful for the discussion of the 
role of non-local exchange, upon which a clean analysis 
of the effects of (dynamical) electronic correlation can be 
then established. 
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